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ABSTRACT 

We present results from multi-epoch spectral analysis of XMM-Newton and Chandra observations of 
the broad absorption line (BAL) quasar APM 08279+5255. Our analysis shows significant X-ray BALs 
in all epochs with rest-frame energies lying in the range of ~ 6.7-18 keV. The X-ray BALs and 0.2-10 keV 
continuum show significant variability on timescales as short as 3.3 days (proper time) implying a source 
size-scale of ~ 10 r g , where r g is the gravitational radius. We find a large gradient in the outflow velocity 
of the X-ray absorbers with projected outflow velocities of up to 0.76 c. The maximum outflow velocity 
constrains the angle between the wind velocity and our line of sight to be less than ~ 22°. Based on our 
spectral analysis we identify the following components of the outflow: (a) Highly ionized X-ray absorbing 
material with an ionization parameter in the range of 2.9< log£<3.9 (the units of £ are erg cm s _1 ) 
and a column density of logiVn ~ 23 (the units of arc cm~ 2 ) outflowing at velocities of up to 
0.76 c. (b) Low-ionization X-ray absorbing gas with logTVu ~ 22.8. We find a possible trend between 
the X-ray photon index and the maximum outflow velocity of the ionized absorber in the sense that 
flatter spectra appear to result in lower outflow velocities. Based on our spectral analysis of observations 
of APM 08279+5255 over a period of 1.2 years (proper time) we estimate the mass-outflow rate and 
efficiency of the outflow to have varied between 16^g 2 M Q yr^ 1 and 64^40 Mq yr _1 and 0.181°'^ to 
l-Ttj. a> respectively. Assuming that the outflow properties of APM 08279+5255 are a common property 
of most quasars at similar redshifts, our results then imply that quasar winds are massive and energetic 
enough to influence significantly the formation of the host galaxy, provide significant metal enrichment to 
the interstellar medium (ISM) and intergalactic medium (IGM), and are a viable mechanism for feedback 
at redshifts near the peak in the number density of galaxy mergers. 

Subject headings: galaxies: active — quasars: absorption lines — quasars: 

individual (APM 08279+5225) — X-rays: galaxies — gravitational lensing 



1. INTRODUCTION 

A number of physical mechanisms of interactions of 
AGNs with their environments have been proposed to ex- 
plain feedback in galaxies and clusters of galaxies. These 
feedback mechanisms include radio jets, AGN winds, and 
AGN radiative heating. It may be the case that several of 
these mechanisms operate simultaneously and their rela- 
tive contribution to the feedback process varies depending 
in part on the size of the super- massive black hole (SMBH), 
and the combined gravitational potential of the SMBH and 
host galaxy. The gas and dust content of the host galaxy, 
the redshift and age of the host, and the inflow and outflow 
properties of the central super-massive black hole may be 
important as well. 

The radio-jet feedback mechanism was proposed after 
the discovery of significant radio cavities and shock fronts 
in several clusters of galaxies. The cavity sizes are large 
enough to imply that the amount of energy injected in the 
intracluster medium (ICM) , as inferred from the mechani- 
cal work pdV required to displace the gas in the radio cav- 
ities, is sufficient to balance radiative losses, suppress cool- 
ing flows, and halt star formation (e.g., Fabian et al. 2009, 
and references therein). The radio jets are thought to be 
driven by either accretion onto the black hole (e.g., Bland- 
ford & Payne 1982) or the black- hole spin (e.g., Blandford 



& Znajek 1977). The details of the process by which radio 
jets deposit energy and distribute it uniformly throughout 
the ICM is currently not well understood. 

Another possibly important mode that may contribute 
to feedback is AGN radiation. Radiative heating (e.g., 
Ciotti & Ostriker 2007) has been proposed as an alterna- 
tive to the radio-jet mechanism to explain the suppression 
of cooling flows in isolated elliptical galaxies and possibly 
in clusters of galaxies. 

In field galaxies, especially ones in the redshift range of 
z « 1 — 3 where the number density of galaxy mergers 
is thought to peak (e.g., Di Matteo et al. 2005), quasar 
winds are thought to be one of the major contributors 
to feedback. Winds can be driven by radiation pressure, 
magneto-centrifugal forces, and thermal pressure or a com- 
bination of these processes. The potential importance of 
quasar outflows has been explicitly demonstrated in the- 
oretical models of structure formation and galaxy merg- 
ers that incorporate the effects of quasar outflows (e.g., 
Silk & Rees 1998; Scannapieco & Oh 2004; Granato et al. 
2004; Springel, Di Matteo, & Hernquist 2005; Hopkins et 
al. 2005, 2006). Observational evidence to support these 
theories of quasar outflows was provided by the discovery 
of near-relativistic X-ray absorbing winds in several broad 
absorption line (BAL) quasars (e.g., Chartas et al. 2002, 
2003) and later in several non-BAL quasars (e.g., Reeves 



1 Department of Physics and Astronomy, College of Charleston, Charleston, SC, 29424, USA, chartasg@cofc.edu 

2 Department of Astronomy & Astrophysics, Pennsylvania State University, University Park, PA 16802, 



chartas@astro.psu.edu, 

saez@astro.psu.edu, niel@astro.psu.edu, giustini@astro.psu.edu, garmire@astro.psu.edu 

3 Dipartimento di Astronomia, Universita degli Studi di Bologna, via Ranzani 1, 40127 Bologna, Italy, margharita.giustini2@iasfbo.inaf.it 

4 INAF/Istituto di Astrofisica Spaziale e Fisica Cosmica di Bologna, via Gobetti 101, 40129 Bologna, Italy, margherita.giustini@unibo.it 



2 



CHARTAS ET AL. 



et al. 2003; Pounds et al. 2003, 2006; Braito et al. 2007). 

The relative contribution of each of these proposed feed- 
back mechanisms will vary in part depending on the fol- 
lowing conditions: 

(a) The importance of the radio-jet mode will depend on 
the radio loudness and the duty cycle of the central AGN. 
The presence of a radio jet in an AGN is likely to depend 
on the available supply of infalling material to feed the 
black hole and/or the spin of the black hole. We note that 
several observations imply that the fraction of radio-loud 
objects varies with redshift and luminosity (e.g., Peacock 
et al. 1986; Schneider et al. 1992; La Franca et al. 1994; 
Jiang et al. 2007). Specifically it is found that the radio- 
loud fraction decreases with redshift and increases with 
luminosity. The radio jet mode therefore may be more 
important at later times and in clusters of galaxies and 
in nearby massive ellipticals that have a larger probability 
of containing a radio-loud object compared to galaxies at 
z > 1 which are more likely to contain a radio-quiet AGN. 

(b) The feedback process is also expected to depend on 
the mass of the central black hole and the gravitational 
potential of the host galaxy. In order for a feedback mech- 
anism to be effective it must be able to drive gas out of the 
host with a velocity that is larger than the escape velocity 
determined by the gravitational potential of the system. 
The mass outflow rates and outflow efficiencies of AGN 
winds will also depend strongly on the mass and luminos- 
ity of the SMBH. Analyses of the kinematic properties of 
UV absorbers in quasars indicate a strong increase in out- 
flow velocity of UV absorbers with quasar UV luminosity 
(e.g., Laor & Brandt 2002; Ganguly & Brotherton 2008). 
Luminous quasars are found to contain winds with veloci- 
ties of the UV outflowing material of up to 60,000 km s _1 
whereas Seyfert galaxies typically have winds with signif- 
icantly lower velocities of up to ~ 2,000 km s _1 . 

(c) Hydrodynamic simulations of AGN atmospheres 
driven by radiation pressure indicate that the strengths 
of the outflows originating from AGN accretion disks de- 
pend on the Eddington ratio (iBol/^Edd) with more mas- 
sive and faster winds produced at larger Eddington ratios 
(e.g., Proga, Stone & Drew 1988; Proga, Stone & Kallman 
2000; Proga & Kallman 2004). In systems with a limited 
supply of gas to fuel the AGN it is expected that AGN 
winds will be relatively weak. 

(d) In the case of a galaxy merger the degree to which 
each mechanism contributes to the overall feedback pro- 
cess may change with time since the merger event oc- 
curred. 

Observations in X-rays of the BAL quasar 
APM 08279+5255, the mini-BAL quasar PG 1115+080, 
and perhaps the low-ionization BAL quasar H 1413-117 
have strongly suggested the presence of near-relativistic 
outflows of X-ray absorbing material in these objects 
(Chartas et al. 2002, 2003, 2007a, 2007b). Observations of 
these quasars in the optical and UV indicate outflow veloc- 
ities and column densities that are more than an order of 
magnitude lower than those inferred from the X-ray spec- 
tra. These differences in outflow velocities and absorbing 
column densities suggest that X-ray BALs probe a highly 
ionized, massive and high- velocity component of the wind 

The WWW sites listing the updates are located at 
http : //xmm2 . esac . esa. int/external/xmm_sw_cal/calib/ respectively. 



that is largely distinct from the absorbers detected in the 
optical and UV wavebands. Most of the mass and en- 
ergy of quasar winds is likely carried by the outflowing 
X-ray absorbers, and observations in the X-ray band are 
therefore crucial for improving our understanding of the 
contribution of quasar winds to feedback in galaxies. 

This paper will focus on determining the significance of 
the wind in the z — 3.91 gravitationally lensed BAL quasar 
APM 08279+5255 in the feedback process during an epoch 
near the peak of the number density of mergers. Specif- 
ically, we present results from analyses of several XMM- 
Newton and Chandra observations of APM 08279+5255 in 
order to determine the outflow velocity and its variability, 
obtain insights on the acceleration mechanism, constrain 
the geometry of the wind, and estimate the mass outflow 
rate and the outflow efficiency. 

Throughout this paper we adopt a A-dominatcd cosmol- 
ogy with H = 70 km s _1 Mpc -1 , il A = 0.7, and fi M = 
0.3. 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. Observations and Data Reduction 

APM 08279+5255 was observed with XMM-Newton 
(Jansen et al. 2001) on 2002 April 28, 2007 October 6, 
and 2007 October 22, respectively. It was also observed 
with the Advanced CCD Imaging Spectrometer (ACIS; 
Garmire et al. 2003) on board the Chandra X-ray Obser- 
vatory (hereafter Chandra) on 2002 February 24 and 2008 
January 14. 

The spectral analysis of the 2002 Chandra observation 
of APM 08279+5255 has been presented in Chartas et al. 
(2002), and the spectral analysis of the 2002 XMM-Newton 
observation of APM 08279+5255 in Hasinger et al. (2002). 
Because of recent significant improvements in the cali- 
bration of the instruments on board Chandra and XMM- 
Newton since the publication of the APM 08279+5255 re- 
sults, we have re-analyzed all observations. Updates on the 
calibration of Chandra and XMM-Newton are reported on 
the Chandra X-ray Center (CXC) and XMM-Newton Sci- 
ence Operations Centre (SOC) World Wide Web (WWW) 
sites, respectively^ 

We analyzed the XMM-Newton data for 
APM 08279+5255 with the standard analysis software 
SAS version 8.0.0 provided by the XMM-Newton SOC. 
The Chandra observations of APM 08279+5255 were an- 
alyzed using the standard software CIAO 4.1 provided by 
the CXC. A log of the observations that includes obser- 
vation dates, observed count rates, total exposure times, 
and observational identification numbers is presented in 
Table 1. 

We note that the count rate and total number of counts 
(~ 40,918 counts from all XMM-Newton observations) for 
APM 08279+5255 are the highest of any BAL quasar X- 
ray spectrum observed to date. This can be compared to 
the typical total number of 0.5-8 keV counts detected in 
BAL quasars obtained from the Sloan digital sky survey 
data release 5 that ranges between a few to ~ 200 counts 
(Gibson et al. 2009). 

For the reduction of the XMM-Newton observations we 
filtered the pn (Striider et al. 2001) and MOS (Turner 

http: //asc .harvard. edu/ciao/releasenotes/history .html and 
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et al. 2001) data by selecting events corresponding to in- 
strument PATTERNS in the 0-4 (single and double pixel 
events) and 0-12 (up to quadruple pixel events) ranges, re- 
spectively. Several moderate-amplitude background flares 
were present during the XMM-Newton observations of 
APM 08279+5255. The pn and MOS data were filtered 
to exclude times when the full-ficld-of-view count rates 
exceeded 20 cnts s _1 and 4 cnts s _1 , respectively. The 
extracted spectra from the pn and MOS were grouped to 
obtain a minimum of 100 counts in each energy bin, al- 
lowing use of x 2 statistics. Background spectra for the 
pn and MOS detectors were extracted from source-free re- 
gions near APM 08279+5255. 

The pn and MOS spectra were fit with a variety of mod- 
els employing XSPEC version 12.5 (Arnaud 1996). The en- 
ergy ranges used for fitting the pn and MOS data were 
0.2-10 keV and 0.4-10 keV, respectively. We performed 
spectral fits to the pn spectra alone and to the pn and MOS 
data simultaneously. Both approaches resulted in values 
for the fitted parameters that were consistent within the 
errors, however, the fits to the higher quality pn data alone 
provided higher quality fits as indicated by the reduced \ 2 
values of these fits. We therefore consider the results from 
the fits to the pn data alone more reliable especially for 
characterizing the properties of the X-ray absorption fea- 
tures. 

For the reduction of the Chandra observations we used 
standard CXC threads to screen the data for status, grade, 
and time intervals of acceptable aspect solution and back- 
ground levels. The pointings placed APM 08279+5255 on 
the back-illuminated S3 chip of ACIS. To improve the spa- 
tial resolution we removed a ± 0"25 randomization applied 
to the event positions in the CXC processing and employed 
a sub-pixel resolution technique developed by Tsunemi et 
al. (2001). 

In both the XMM-Newton and Chandra analyses we 
tested the sensitivity of our results to the selected back- 
ground and source-extraction regions by varying the loca- 
tions of the background regions and varying the sizes of 
the source-extraction regions. We did not find any signifi- 
cant change in the background-subtracted spectra. For all 
models of APM 08279+5255 we included Galactic absorp- 
tion due to neutral gas with a column density of JVh=3.9 
x 10 20 cm~ 2 (Stark et al. 1992). All quoted errors are at 
the 90% confidence level unless mentioned otherwise. 

2.2. Chandra and XMM-Newton Spectral Analysis of 
APM 08279+5255 

We first fitted the Chandra and XMM-Newton spectra 
of APM 08279+5255 with a simple model consisting of a 
power law with neutral intrinsic absorption at z = 3.91 
(model 1 of Table 2). These fits are not acceptable in a 
statistical sense as indicated by the reduced x 2 ■ The resid- 
uals between the fitted simple absorbed power-law model 
and the data show significant absorption for energies in 
the observed-frame band of < 0.6 keV (referred to hence- 
forth as low-energy absorption) and 2-5 keV (referred to 
henceforth as high-energy absorption). 

To illustrate the presence of these low- and high-energy 
absorption features, we fit the spectra from observed-frame 
4.5-10 keV with a power-law model (modified by Galac- 

6 The equivalent width (EW) is defined as EW = J Fc ~ Fe dE, where 



tic absorption) and extrapolated this model to the energy 
ranges not fit. The residuals of these fits are shown in Fig- 
ures 1 and 2. Significant low- and high-energy absorption 
is evident in all observations. 

We proceed by fitting the spectra of APM 08279+5255 
with the following models: 1) absorbed power-law (APL), 
2) absorbed power-law with a notch (APL+No), 3) 
ionized- absorbed power-law with a notch (IAPL+No), 4) 
absorbed power-law with an absorption edge (APL+Ed), 
5) ionized-absorbed power-law with an absorption edge 
(IAPL+Ed), 6) absorbed power-law with two absorption 
lines (APL+2AL), 7) ionized-absorbed power-law with two 
absorption lines (IAPL+2AL), 8) absorbed power-law with 
two intrinsic ionized absorbers (APL + 2IA), and 9) ab- 
sorbed power-law with with two partially covered intrinsic 
ionized absorbers (APL+PC*(2IA)). The XSPEC nota- 
tions for these models are given in the notes of Tables 2 
and 3. 

The results from fitting these models to the three XMM- 
Newton and two Chandra spectra are presented in Tables 
2 and 3. For spectral fits using models 3, 5, and 7 the low- 
energy absorption is modeled using the photo-ionization 
model absori contained in XSPEC (Done et al. 1992). We 
note that the absori model is just a first approximation 
to what is likely a more complex situation. As an inde- 
pendent check of the accuracy of the fits that used the 
absori model we also repeated several of these fits using 
the warm-absorber model XSTAR (Kallman et al. 1996; 
Kallman & Bautista 2001). XSTAR calculates the physical 
conditions and emission spectra of photoionized gases. In 
the current analysis we use a recent implementation of the 
XSTAR model that can be used within XSPEC. 

In most epochs, models APL+2AL and IAPL+2AL 
that included two absorption lines provided significant im- 
provements to the fits compared to models that included 
an absorption edge. Specifically, the F-test indicates 
the significance of the improvements of the IAPL+2AL 
model compared to the IAPL+ Ed model to be > 99.99% 
(A X 2 = 42.5), >84% (A X 2 = 6.4), >74% (A x 2 = 6.3), 
>99% (Ax 2 = 15.6), and >99.7% (A X 2 = 17.1), for four 
additional parameters in epochs 1 through 5, respectively. 

We note that the X-ray BALs in the 1.5-3.6 keV band 
(observed-frame) correspond to a significant detection fol- 
lowing the criteria described in §3 of Vaughan & Uttley 
(2008). Specifically, we find the ratio of the total equiva- 
lent widths^ of the absorption features to their uncertainty 
to be EW/ctew ^ 3 in every observation (see models 6 and 
7 in Table 2). 

As indicated in the results of the spectral fits shown in 
Table 2, the energies and widths of the absorption lines of 
models APL+2AL and IAPL+2AL vary significantly be- 
tween several epochs. This is also suggested by the resid- 
uals shown in Figures 1 and 2, where we indicate with 
arrows the locations of several of the best-fit energies of 
the absorption lines. 

To estimate the significances of variations of the best-fit 
energies (-E a bs) and widths (cr a bs) of the absorption lines 
between epochs we calculated the \ 2 confidence contours 
of these best-fit parameters. 

We first focus on epochs 1 and 5 where the high-energy 
absorption lines appear to be narrower than in other 

F c is the continuum flux and Fe is the flux in the absorber. 
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epochs (see Figures 1 and 2) and also show variability in 
energy and width. In Figure 3 we show the 68%, 90%, 
and 99% \ 2 confidence contours of -E a b s versus normaliza- 
tion of the first and second absorption lines (APL + 2 AL 
model) in epoch 1 (top panel) and epoch 5 (bottom panel). 
Since the 99% and 95% confidence contours of the second 
absorption line (component abs2) between epochs 1 and 
5, respectively, do not overlap in energy we conclude that 
a significant change occurred at the > 99.9% confidence 
level of the energy of the second absorption line between 
these epochs. 

In Figure 4 we show the 68%, 90% and 99% \ 2 confi- 
dence contours of versus cr a b s of the first and second 
absorption lines (APL + 2 AL model) in epochs 1 and 5. 
The probability that the energy-width parameters of the 
first absorption line are the same between epochs 1 and 5 
(null probability) is <1 x 10" 4 []. 

An additional significant variation of particular interest 
is found between epochs 3 and 4 that are separated by 
only 3.3 days in the rest frame. The 0.2-10 keV pn count 
rate of APM 08279+5255 varied between epochs 3 and 4 
by 36.8 ± 0.3 %. We investigated whether the X-ray BALs 
varied in this 3.3 day period. In Figure 5 we show the 68%, 
90% and 95% x 2 confidence contours of E a b s versus cr a b s of 
the first absorption line (APL + 2 AL model) in epoch 3 
(dotted line) and epoch 4 (solid line). The 95% confidence 
contours of epoch 3 almost touch the 90% confidence con- 
tours of epoch 4. The probability that the energy-width 
parameters of the first absorption line (model 7; Table 2) 
are the same between epochs 3 and 4 (null probability) is 
< 0.05x0.10 = 5 x 10~ 3 . 

In Table 3 we present results from spectral fits that at- 
tempt to model the high-energy absorption with two out- 
flowing ionized absorbers that are either completely cover- 
ing the source (model 8) or are partially covering it (model 
9). For spectral fits using models 8 and 9 the high-energy 
absorption in modeled using XSTAR. The low-energy ab- 
sorption is modeled with a stationary neutral absorber. 
For epochs 1 and 5 we find that the combined covering fac- 
tor of the highly ionized absorbers is close to unity. This 
is expected since the detected ionized absorption lines in 
epochs 1 and 5 are narrower than the ones detected in the 
other epochs whereas the presence of any significant di- 
rect component (emission that does not get absorbed by 
the outflowing material) would have diluted any narrow 
absorption lines. For epochs 2, 3 and 4 the spectral fits 
(model 9 of Table 3) do allow for an emitted component 
that is not absorbed by the ionized outflowing material. 

Our fits with XSTAR indicate an outflowing X-ray ab- 
sorbing medium with ionization parameters ranging be- 
tween log£ = 2.9^1 and log£ = 3.9 (model 8 of Table 3). 
The two strongest iron lines for this highly ionized absorb- 
ing medium have rest (or laboratory) energies of 6.70 keV 
(Fe xxv Is 2 - \s2p) and 6.97 keV (Fe xxvi Is - 2p). In 
general the Fe xxv Is 2 — ls2p line will be stronger than the 
Fe xxvi Is - 2p line for a medium with 2.75 < log £< 4.0 



(see Figure 3 of Saez et al. 2009). 

We emphasize that the photo-ionization models used in 
our analysis (i.e., absori and XSTAR) do not consider possi- 
ble velocity gradients in the outflowing absorber and there- 
fore cannot provide realistic models of the X-ray BALs. 
We attempt to mimic the velocity broadening of the lines 
by introducing in the XSTAR model large turbulent veloci- 
ties of Wtmb = 10,000 km s _1 for epochs 1 and 5 and iWb 
= 30,000 km s -1 for epochs 3, 4 and 5. 

3. DISCUSSION 

3.1. Short- and Long-Term Variability of 
APM 08279+ 5255 

As we showed in §2 the 0.2-10 keV pn count-rate of 
APM 08279+5255 varied by - 36.8 ± 0.3 % between 
epochs 3 and 4. From a light-travel time argument we 
estimate that the observed short-term variability between 
epochs 3 and 4 implies a size-scale of the X-ray emis- 
sion region of the order of / V ar = cAt/(l + z) ~ 7.4 x 
10 15 cm. We compare this emission-region size scale to 
the radius of the innermost stable circular orbit (nsco) 
which for Schwarzschild and Kerr (maximally spinning) 
black holes are 6r g and r s , respectively, where r g = GM/c 
is the gravitational radius. For APM 08279+5255 we as- 
sume Mbh ~ 10 12 M©/i£ , where \xl ~ 100 (Egami et al. 
2000)0 We find that nsco = 4.5 x 10 15 cm for the case 
of a Schwarzschild black hole which is comparable to the 
size-scale of the X-ray emitting region implied by our light- 
travel time argument. This result is consistent with recent 
X-ray and UV (rest-frame) monitoring observations of the 
z — 2.32 gravitationally lensed quasar HE 1104—1805 that 
have constrained the size of the X-ray emitting region in 
this object to be smaller than 6r g and the size of the UV 
emission region to be ~ 30r g based on the analysis of the 
microlensed light-curves of the lensed images (Chartas et 
al. 2009). A hint to the cause of the spectral variability 
can be seen in Figure 6 where we show the three over- 
plotted XMM-Newton spectra of APM 08279+5255. 

The spectra between epochs 3 and 4 appear to differ 
over the entire observed energy band by approximately a 
constant. This type of behavior suggests that the cause of 
the variability is not the result of a change in the absorb- 
ing column density but is likely due to either a change in 
covering factor or a change in the ionizing flux from the 
central source. 

As shown in Figure 6 the comparison between the epoch 
2 and 3 XMM-Newton spectra, that are separated by ~ 
5.4 years (observed- frame) , indicates that they converge 
at energies above ~ 3 keV but diverge significantly be- 
low this value. This type of energy-dependent variation 
is indicative of a decrease in the absorbing column den- 
sity between epochs 2 and 3 and is consistent with the 
results of our spectral analysis shown in Table 2. Since 
the relevant energy range (> 15 keV in the rest frame) 
is not significantly affected by neutral or ionized absorp- 
tion, we interpret the variation between the epoch 2 and 3 

7 The product of the probabilities of being outside the confidence contours that barely touch (see Figure 4) is an upper limit to the null 
probability 

8 Throughout this paper we adopt the definition of the ionization parameter of Tarter et al. (1969) given by £ = ^"^ = /jj 000 ""'' 
where nn is the hydrogen number density, and r is the source-cloud separation. 

9 See, however, Riechers et al. (2009), who find a magnification of fXL 
to obtain a black-hole mass of Mbh ~ lO 11 /^ Mq. 
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4. Riechers et al. (2009) also use the observed width of the CIV line 
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XMM-Newton spectra as an indication that the intrinsic 
unabsorbed luminosity of APM 08279+5255 has not sig- 
nificantly varied between these two observations and that 
most of the observed changes in flux result from changes in 
the opacities of the low and/or high-ionization absorbers. 

In Figure 7 we over-plot the epoch 1 and epoch 5 Chan- 
dra spectra that are separated by <~ 5.9 years (observed- 
frame). The spectra appear to converge at energies be- 
low ~ 0.7 keV and above ~ 5 keV and significantly differ 
between these energies. As shown in our spectral anal- 
ysis, the X-ray BALs varied significantly between these 
two epochs which is consistent with the observed differ- 
ence of the spectra between 0.7-3.6 keV seen in Figure 7. 
The variability of the absorption between epochs 1 and 
5 can in general be caused by at least three factors: (a) 
transverse motion of the absorber across the continuum 
source, (b) changes in the opacity of the absorber, and (c) 
changes in the velocity structure of the absorber. Epochs 

1 and 5 are separated by 2150 days which is significantly 
longer than the expected timescale of variability due to 
transverse motion and also longer than the expected flow 
timescale. We therefore cannot reach a firm conclusion re- 
garding the cause of the variability between epochs 1 and 
5. The significant differences in timescales, however, imply 
that it is unlikely that the variability is being produced by 
the same absorber moving across our line of sight. 

To better illustrate the flux variability of 
APM 08279+5255 we show the ratios of the observed 
and best-fit modeled (using model 6 of Table 2) flux den- 
sities Fj /Fj between the three XMM-Newton observations 
in Figure 8 and between the two Chandra observations in 
Figure 9, where i and j represent the epochs compared, 
F is the flux density in units of counts s _1 keV -1 and 
At is the observed-frame time between epochs. We first 
focus on the short-term variability between epochs 3 and 
4 shown in Figure 8. We notice that the variability of 
the flux density is to first order independent of energy 
over the observed 0.5-10 keV band with the exception of 
apparent variability in the region of 1.5-3.6 keV (observed- 
frame) near the Fe BALs. A plausible explanation of the 
energy independence of the short-term variability of the 
flux density may involve motion of the highly ionized ab- 
sorbing outflow perpendicular to our line of sight. This 
motion may change the covering factor resulting in possi- 
ble energy-independent variability of the flux density. A 
second possible explanation is a change in the ionizing 
flux of the central source. Our spectra fits that assumed 
a model that includes emission through two ionized ab- 
sorbers and emission that is not absorbed by the ionized 
absorber (model 9 of Table 3) suggest a possible decrease 
in the total covering factor of the ionized absorbers be- 
tween epochs 3 and 4. Unfortunately the limited S/N of 
the current spectra provide poor constraints on the cov- 
ering factors and therefore the implied change in covering 
factor between epochs 3 and 4 can only be considered 
suggestive. 

We second center on the long-term variability between 
epochs 2 and 3 shown in Figure 8. The variability of the 
flux density between these epochs becomes large below ~ 

2 keV, however, no significant variability is detected above 
this energy. A possible explanation of such energy depen- 
dent variability is a decrease in column density between 



epochs 2 and 3. 

We finally comment on the long-term variability be- 
tween epochs 1 and 5 shown in Figure 9. We note that 
the flux density appears to have varied by a larger fac- 
tor near the region of the Fe BALs that extends above ~ 
1.5 keV (observed frame). A possible explanation for this 
flux-density variability between these epochs is a change 
in the properties of the outflowing ionized absorbers. This 
interpretation is consistent with the results of our previous 
spectral analysis that showed a significant change in the 
energies and widths of the absorptions lines (see Figures 3 
and 4) . 

3.2. Mass-Outflow Rate and Efficiency of the Outflow 

Our XMM-Newton and Chandra observations of 
APM 08279+5255 provide constraints on the velocity, col- 
umn density, and location of the X-ray BAL material that 
allow us to estimate the mass-outflow rate from the ex- 
pression: 

Mi = 4:nri(ri/Ari)N Ki impV wiru i,ifc,i- (1) 

where Ar; is the thickness of the absorber at radius r;, iVn,i 
is the hydrogen column density, w w i n d,i is the outflow veloc- 
ity of the X-ray absorber, / Cj ; is the global covering factor 
of the absorber, and i indicates the absorbing component. 
Combining the constraints on the mass-outflow rate and 
the bolomctric luminosity we can estimate the efficiency of 
the quasar outflow. The efficiency is defined as the ratio 
of the rate of kinetic energy injected into the interstellar 
medium (ISM) and intergalactic medium (IGM) by the 
outflow to the quasar's bolomctric luminosity, i.e., 

1 M-iP ifi 

1 i^wind,i o -f ( /A \ t\j wind,i ln \ 

£k,i = — ; = 27r/ Ci iri(ri/Ari)iVH,iTOp— (2) 

^ -^Bol ^Bol 

where LboI is the bolometric luminosity of the quasar. 

The best-fit values of the equivalent widths of the X-ray 
absorption lines (assuming model 6 of Table 2) were used 
to estimate the hydrogen column densities JVh of the X- 
ray BALs using a curve-of-growth analysis (e.g., Spitzer 
1978). For our crude calculation we assumed b parame- 
ters of the order of the observed widths of the lines where, 
b = \/2er u and a u is the velocity width of the line. The bulk 
outflow velocities of each outflow component were inferred 
from the best-fit values of the energies of the blucshifted 
absorption lines obtained by spectral fits to the data as- 
suming model 6 of Table 2. To convert from line energies 
to outflow velocities we used the relativistic Doppler for- 
mula assuming that the high-energy absorption lines are 
the result of resonance absorption from ions of Fe xxv in 
a gas with solar abundances, and the angle between our 
line of sight and the outflow direction is zero. For estimat- 
ing the outflow efficiency of APM 08279+5255 we assumed 
Lboi = 7 x 10 15 ^£ 1 L Q (Irwin ct al. 1998; Reichers et al. 
2009). We assume a conservatively wide range for the cov- 
ering factor of f c =0. 1-0.3 based on the observed fraction 
of BAL quasars (e.g., Hewett & Foltz 2003). We obtained 
a constraint on the covering fraction of the flow by plac- 
ing limits on the presence of any emission components of 
the lines. Crudely the ratio of the emission to absorp- 
tion equivalent widths is proportional to the global cover- 
ing factor f c . The tightest upper bound is found for the 
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XMM-Newton observation of APM 08279+5255 in epoch 
3. We find that EW eml / EW ahsl < 0.24 at the 90% con- 
fidence level assuming model 6 of Table 2 and we expect 
f c oc EWeml /-EWabsi- We note that the covering factor 
f c will depend on the fluorescence yields of the ionized 
absorber which for Fe xxv and Fe xxvi are 0.5 and 0.7, 
respectively ( e.g., Krolik & Kallman 1987). For this crude 
estimate we included a narrow emission line in model 6 at a 
rest-frame energy of 6.7 keV. The exact relevant energy of 
the emission line to consider is probably geometry depen- 
dent. For estimating the outflow efficiency we assumed 
a fraction r/Ar ranging from 1-10 based on theoretical 
models of quasar outflows (e.g., Proga & Kallman 2000). 
Assuming that the maximum outflow velocity is produced 
by gas that has reached its terminal velocity one obtains 
the approximation i?i aunc h ~ -R s (c/w w ind) 2 , where v wind 
is the observed outflow velocity. Based on our estimated 
maximum outflow velocity (v = 0.76 c) we expect r to 
be similar to i?i aun ch arid range between 3 R s and 15 R s 
(see Figure 4 of Saez et al. 2009). We note that the short 
term variability time-scale of the X-ray BALs (3.3 d in 
rest-frame) is consistent with a launching radius of a few 
times Rs. 

We used a Monte Carlo approach to estimate the errors 
of M; and ek- The values of / c , r/Ar, and n were assumed 
to have uniform distributions within their error limits. The 
values of Nh were assumed to have normal distributions 
within their error limits. By multiplying these distribu- 
tions and with the appropriate constants from equations 1 
and 2 we obtained the distributions of M and et- We fi- 
nally determined the means of the distributions of Mi and 
ek and estimated the 68% confidence ranges. 

In Table 4 we present the outflow velocities, the column 
densities, the mass-outflow rates and the efficiencies of the 
outflowing components. 

3.3. Driving Mechanism of X-ray BALs 

One of the major unanswered questions in current the- 
oretical and numerical models of quasar winds involves 
explaining how highly ionized outflows of X-ray absorbing 
material become accelerated to near-relativistic velocities. 
The main problem is that for radiatively driven winds the 
magnitude of the force on an absorber is a function of its 
ionization state. As the ionization parameter of the ab- 
sorber increases, fewer resonant transitions are available to 
absorb photons from the source, thus resulting in a weaker 
driving force. Chelouche & Netzer (2001) have explored 
the dependence of the radiation driving force on the ion- 
ization parameter of the outflowing gas and they conclude 
that for a highly ionized X-ray absorbing gas, the average 
force multipliers are close to 10. 

Most current theoretical work on radiative acceleration 
in quasars has focused on interpreting the outflows of UV 
absorbers (e.g, Murray et al. 1995; Proga et al. 2000; 
Proga et al. 2004; Everett 2005), however; there are no 
self-consistent models that can explain the acceleration 
mechanism that leads to the near-relativistic outflows of 
X-ray absorbing material. 

A clue to understanding the acceleration process of the 
wind in APM 08279+5255 is perhaps provided by the ob- 
served large difference between the maximum velocities of 

10 The dimensionless optical depth is t = n e cTTV t ^/ (dv/dr), where, n e 



the UV (vuv ~ 0.04c; e.g., Srianand & Petitjean 2000) 
and X-ray (vx-ray ~ 0.43-0. 76c) absorbers. One possible 
explanation of the difference between the UV and X-ray 
outflow velocities is that the UV and X-ray BALs are pro- 
duced by different absorbers, with UV emission from the 
accretion disk driving the UV absorbers and X-ray emis- 
sion from the hot corona contributing to the acceleration 
of the X-ray absorbers. If the X-ray absorbers are partly 
driven by radiation from the hot corona we might expect 
to detect a correlation between the properties of the X-ray 
BALs and the properties of the X-ray spectrum. For ex- 
ample, we might expect the maximum outflow velocity of 
the X-ray absorbers to depend on the shape of the X-ray 
spectrum that at first order can be represented with the 
X-ray photon index (r). 

The minimum and maximum projected velocities (v m i n , 
t) mm ) of the outflow are estimated from the minimum 
and maximum energy ranges (-E mm , i? max ) of the high- 
energy absorption features in APM 08279+5255 assum- 
ing that the high-energy X-ray BALs are produced by the 
resonance transition Fe xxv (Is 2 — ls2p). We obtained 
-E m in and £ max from our spectral fits assuming the two 
absorption-line (APL+2AL) model. Specifically, based on 
the best- fit values of the APL+2AL model (model 6; Ta- 
ble 2), we obtain E min = E ahsl - 2cr abs i and E max = 
E a hs2 + 2a a bs2 • The values of E min and E max are presented 
in Table 5. 

In Figure 10 we show the maximum outflow velocity 
observed from APM 08279+5255 as a function of X-ray 
photon index based on the XMM-Newton, Chandra, and 
Suzaku observations. We notice a possible trend between 
the photon index and the maximum outflow velocity in the 
sense that flatter X-ray spectra appear to result in lower 
outflow velocities. One possible explanation is that flatter 
X-ray spectra over-ionize the X-ray absorber resulting in 
a decrease of the force multiplier and a lower outflow ve- 
locity. If this trend is confirmed with future observations 
it would imply that the dominant acceleration mechanism 
responsible for the near-relativistic velocities of the X-ray 
outflowing gas in BAL quasars is indeed radiation driving 
as opposed to magnetic driving. 

We investigated our hypothesis of the origin of the possi- 
ble T — w max correlation by calculating the force multiplier 
as a function of the incident spectral energy distribution 
(SED). The force multiplier represents the ratio by which 
the bound-bound (line) and bound-free (continuum) opac- 
ity increases the radiation force relative to that produced 
by Thomson scattering alone. 

We performed calculations of force multipliers assuming 
a thin slab illuminated by an ionizing continuum using the 
CLOUDY code. We compared our results to those of Arav 
et al. (1994) and Everett et al. (2005) who have used a 
similar approach to ours. In general our computations are 
in good agreement. We performed calculations assuming 
the SED is a power-law extending from the UV (or 1 Ryd 
-13.6 eV) to hard X-rays (or 10 4 Ryd -100 keV). For 
our calculations we assume SEDs with two different val- 
ues of the power-law photon index. A "soft" (r = 2.1) 
and a "hard" (r =1.7) SED. For each SED we calculated 
the continuum (Mq) and the line (Ml) components of the 
force multiplier. Ml depends on an additional parame- 

is the electron number density, <tt is the Thomson cross section and u t h 
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ter, t, which is commonly referred to as the "effective 
electron optical depth" and encodes the dynamical infor- 
mation of the wind in the radiative acceleration calculation 
(see equation 2.5 of Arav & Li 1994). For our calculations 
we have assumed log t = — 7. In Figure 11 we show the 
Mc and Ml components of the force multiplier as a func- 
tion of the ionization parameter and for the soft and hard 
SED cases. Our simple ionized absorber model indicates 
that an increase of the photon-index from 1.7 to 2.1 of 
the incident spectrum will result in a large increase of the 
force multiplier. This result is consistent with the possible 
trend between T and u max shown in Figure 10. 

Theoretical models of AGN outflows (e.g., Murray et 
al. 1995; Chelouche & Netzer 2003) and our observations 
of APM 08279+5255 suggest that the UV BAL gas could 
be shielded from the driving radiation by an absorbing 
medium surrounding the central source. We included the 
effect of such a shield on our calculated values of the force 
multiplier by assuming that the soft and hard SEDs are at- 
tenuated by an absorber with log Ah — 22.8 and the tem- 
perature at the illuminated face of the absorbing shield 
has the value log T[K] = 5.0 (log £ ~ 1.4). As shown 
in Figure 11, Mq decreases in the absorbed-SED case for 
log £ < 2, whereas, Ml increases in the absorbed-SED 
case. We find that Ml increases by almost an order of 
magnitude in the absorbed-SED case when the photon in- 
dex increases from 1.7 to 2.1. In Figure 11 we show that 
in addition to the SED and t parameter, the force mul- 
tiplier depends strongly on the ionization parameter. In 
Figure lib the force multiplier for a shielded outflow drops 
from about 1000 to 500 between log £ of 2-2.8 and drops 
from about 500 to 1 between log £ of 2.8-3, i.e., simi- 
lar to the observed range of the ionization parameter of 
APM 08279+5255. We expect that we observe the ionized 
X-ray absorber of APM 08279+5255 while it has obtained 
its terminal velocity. The observed large value of the ion- 
ization parameter of APM 08279+5255 may therefore not 
be representative of the ionization parameter during the 
initial acceleration phase of the absorber. 

We note that our assumption that the driving force on 
the high-energy absorbers is produced primarily by X-rays 
is plausible since the short term variability time-scale of 
the X-ray BALs of APM 08279+5255 suggest a launching 
radius of < 10i?s (e.g, Chartas et al. 2002; Saez et al. 
2009; this work) and recent studies of AGN employing the 
microlensing technique indicate that the X-ray emission 
region of the hot corona in AGNs is compact with a half- 
light radius of a few i?s and their UV regions are roughly 
a factor of ten larger (e.g, Morgan et al 2008; Chartas 
et al. 2009). Therefore UV radiation is not expected to 
contribute initially at small radii to driving the X-ray ab- 
sorbing outflow. 

If the r — w max trend is confirmed with additional ob- 
servations we plan to produce a more sophisticated model 
that will include more realistic kinematic, ionization and 
absorption properties of the outflow. 

3.4. Constraints on Outflow Direction 

The observed maximum projected outflow velocities of 
the X-ray BALs combined with the constraint that the 
maximum velocity of the X-ray absorber along the outflow 



direction must satisfy w ou tflow < c leads to a constraint on 
the maximum angle between the outflow direction and the 
observed line of sight through the absorber (0 max ). The 
best-fitted values of the ionization parameter for the out- 
flowing absorbers (see model 8 of Table 3) lie in the range 
of 2.9 < log£ < 3.9 and imply that the high-energy absorp- 
tion is the result of resonance absorption from ions of 
Fe xxv and/or Fe xxvi. Since the present X-ray spectra 
cannot distinguish between these ionization stages we esti- 
mate the angle between the outflow direction and our line 
of sight considering that the X-ray BALs are produced by 
resonance transitions of either Fe xxv (Is 2 — ls2p), Fe xxv 
(Is 2 — ls3p), Fe xxvi (Is — 2p) and Fe xxvi (Is - 3p). In 
Figure 12 we show the outflow velocity versus the angle 
between the direction of the outflowing absorber and our 
line of sight through the absorber. The outflow velocities 
were calculated using the relativistic Doppler formula for 
the case of the maximum rest-frame energy of the X-ray 
BALs of E max = 17.9 keV (see Table 5) observed in epoch 
3. For the case where the X-ray BALs are produced by the 
resonance transition of Fe xxv (Is 2 — ls2p) the detected 
projected maximum outflow velocity of v max ~ 0.76c leads 
to the constraint of # max < 22°. Such a small angle is con- 
sistent with the unification scheme of BAL and non-BAL 
quasars that posits that BAL quasars are viewed almost 
along the outflow direction. 

4. CONCLUSIONS 

We have presented results from an analysis of three 
XMM-Newton and two Chandra observations of the z = 
3.91 BAL quasar APM 08279+5255. The main goals 
of these observations were to study the kinematic and 
photoionization properties of the wind in order to assess 
whether it plays an important role in controlling the evo- 
lution of the host galaxy and central black hole. Addi- 
tional objectives included understanding the mechanism 
that drives the X-ray absorbers to near-relativistic veloci- 
ties and constraining the geometry of the outflow. 

The main conclusions of our spectral and timing analy- 
ses are the following: 

(a) X-ray BALs lying in the range of 1.5-3.6 keV (ob- 
served frame) are detected at the greater than 99.9% con- 
fidence level in all five observations of APM 08279+5255. 
We note that a recent analysis of three Suzaku observa- 
tions of APM 08279+5255 also shows significant detections 
of X-ray BALs in all three observations (Saez et al. 2009). 

(b) Based on our spectral analysis we identify the fol- 
lowing components of the outflow. First, a highly ion- 
ized X-ray absorbing material with an ionization param- 
eter in the range of 2.9 < log £< 3.9 and a column density 
of log Ah ~ 23 outflowing at velocities of up to 0.76 c, 
and second a low-ionization X-ray absorbing gas with a 
column density of log Ah ~ 22.8 that may constitutes the 
X-ray shielding gas that is thought to protect the UV ab- 
sorber from becoming over-ionized. Models that include 
partial covering result in even larger column densities of 
the ionized absorber of up to log Ah ~ 24, however, we 
caution that these values are not well constrained (see 
model 9 of Table 3) by the spectral fits. We note that 
no iron overabundance is required to fit the X-ray spectra 
of APM 08279+5255 consistent with the spectral analysis 



is the thermal velocity of the gas. The line force multiplier increases with decreasing t. 
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of previous observations of this object. 

(c) Significant variability of the X-ray continuum and 
X-ray BALs over short and long time-scales is detected. 
Specifically, we detect a ~ 36.8 ± 0.3 % change of the 
0.2-10 keV pn count-rate between epochs 3 and 4 that are 
separated by only 3.3 days (proper-time). Based on a sim- 
ple light-travel time argument this variability time-scale 
implies a radius of the emission region of <~ 7.4 x 10 15 cm 
which is comparable to risco = 4.5 x 10 15 cm for the 
case of a Schwarzschild black hole in APM 08279+5255. 
We speculate based on the energy independent spectral 
change between epochs 3 and 4 (with the exception of the 
Fc B AL region) that the cause of this variability is a change 
in the covering fraction as the outflowing absorbing mate- 
rial moves across our line of sight. Significant changes in 
the flux densities between epochs 2 and 3 that are sepa- 
rated by ~ 5.4 years (observed-frame) are detected. The 
variation of flux density is significant below 3 keV (ob- 
served frame) showing a increase toward lower energies but 
is near zero above 3 keV. This type of energy-dependent 
change suggests that the intrinsic unabsorbed luminosity 
of APM 08279+5255 has not significantly varied between 
these epochs and the long-term variability is the result of 
a decrease in the opacity of the absorber. This interpreta- 
tion is consistent with the results of our spectral analysis. 

(d) Assuming our interpretation that the high-energy X- 
ray BALs are produced by resonance transitions of highly- 
ionized iron we estimate the mass-outflow rate and effi- 
ciency of the outflow to have varied over a period of 1.2 
years (proper-time) between 16^g 2 M Q yr" 1 and 64^q 
M q yr _1 and 0.181°'^ to 1.7±{%, respectively. 

(e) The detected projected maximum outflow velocity 
of fmax ~ 0.76c leads to a constraint on the angle between 
the outflow direction and the observed line of sight through 
the absorber of # max < 22° (assuming Fe xxv(ls 2 — ls2p)) 
and # m ax ~ 26° (assuming Fe xxvi( Is — 2p)). Such a small 
angle is consistent with the unification scheme of BAL and 
non-BAL quasars that posits that BAL quasars are viewed 
almost along the outflow direction. 

(f) A possible trend is found between the X-ray photon 
index and the maximum outflow velocity of the ionized X- 
ray absorber in the sense that flatter X-ray spectra appear 
to result in lower outflow velocities. One possible expla- 
nation is that flatter X-ray spectra over-ionize the X-ray 
absorber resulting in a decrease of the force multiplier of 
the X-ray absorber and thus a lower outflow velocity. 
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Table 1 

Log of Observations of APM 08279+5255 



Observation Date 


Observatory 


Observation 
ID 


Time a 
(ks) 


N b 
net counts 


/o.2-2 C 

erg s" 1 cm -2 


erg s _1 cm" 2 


2002 February 24 (Epoch 1) 
2002 April 28 (Epoch 2) 
2007 October 06 (Epoch 3) 

2007 October 22 (Epoch 4) 

2008 January 14 (Epoch 5) 


Chandra 
XMM-Newton 
XMM-Newton 
XMM-Newton 

Chandra 


2979 
0092800201 
0502220201 
0502220301 

7684 


88.82 
83.46 
56.38 
60.37 
88.06 


5,627 ± 75 
12,820 ± 139 
11,400 ± 114 
16,698 ± 133 

6,938 ± 83 


J-O-0.1 

1 q+oi 

2 5+ 01 

r+0.1 
°-°-0.1 

1 q+ - 2 

L - i >-0.2 


4 3+0.1 
^•°-o.i 

4.i±S:i 
3.9tg:i 

5-0±g;i 



"Time is the effective exposure time remaining after the application of good time-interval (GTI) tables and the removal of portions of the 
observation that were severely contaminated by background flaring. 

'Background-subtracted source counts including events with energies within the 0.2-10 keV band. The source counts and effective exposure 
times for the XMM-Newton observations refer to those obtained with the EPIC PN instrument. See §2 for details on source and background 
extraction regions used for measuring N BC . 

c The absorbed fluxes (in units of 10 -13 ergs cm" 2 s" 1 ) in the 0.2-2 keV and 2-10 keV observed-frame band are obtained using the model 
APL+2AL (model 6; §3). The errors are at the 68% confidence level. 
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Table 2 

Results from spectral fits to the Chandra and XMM-Newton observations of APM 08279+5255. 
Model" Parameter 2 Values Epoch l d Values Epoch 2 d Values Epoch 3 d Values Epoch 4 d Values Epoch 5 d 



r 

iVu 

.2 , 



x7" 



1 74+0.O6 

4 2 o+0.78 

189.8/109 
2.6 x 10~ 6 



1 Q4+0-05 

c Q+0.52 
°- a -0.49 

192.8/120 
0.025 



o ii +0.05 
z - 1± -0.05 



5.26 



+0.40 



191.1/110 

2.8 x 10~ 5 



2 14+ - 04 
5 50 +0 - 31 
244.8/148 
2.6 x 10~ 6 



100.9/75 
9.7 x 10~ 7 



r 

-Enotch [keV] 
W notch [keV] 
fnotch [keV] 
X> 



i 7 r + 0.02 
L - 1 J -0.02 

8.041°;^? 

n 90+O.O6 
u - zo -0 .06 



151/107 
3.2 x 10~ 3 



I. 87 
6.19 

II. 24 



+0.04 
0.05 
+0.51 
0.45 
-0.50 
-0.72 
7 7n +0.90 
' • ' u -1.50 
91 + 04 

122.9/117 
0.34 



2.04Z 
5.391 



,1+0.04 
-0.04 
-,+0.39 
-0.38 

11 62+ ' 44 
ll.OZ_ 3g 

ti.i$f_ 79 
21 + 04 
117.4/107 
0.23 



2.081 
5.641 

12.911 



j+0.04 
-0.03 
4+0.29 
-0.31 
I +1.30 
L -0.91 

n.85t?:» 
o.i6tg:SS 

175.8/145 
0.04 



1 Q5+ 004 

l.»0_g Q 4 



10.85 
8.49 



1.30 
0.96 
+0.2 
0.2 

16+ 007 

1 

83.2/73 
0.19 



r 

iVu 

Enotch [keV] 
W notch [keV] 

fnotch [keV] 

^(x 2 H 



+0.06 
-0.03 
+8.3 
8.6 



1.78 

15.6 
<0.18 
8.0 
0.25 



+0.1 
0.1 
+0.06 
0.06 
1 



146/106 
6.1 x 10~ 3 



1.92: 

13.5j 



5+0.05 
-0.05 
-+4.85 
-4.49 

150±1|? 

ii.i5tg:?? 
8-i2±i;S 

n 1O+0.04 

112.7/116 
0.57 



2.0818:°! 



9.92^ 



,+4.0 

n-54±8:|g 

8 94+ 1 ' 17 

°- ay: -0.73 

20+ 04 

u - zu -0.04 

114.1/106 
0.28 



2 in +0 04 

^•1U_ 04 
,+2.18 
-3.45 



8.36! 



09+97 

12.82l°:S 
11.78li'^ 
l6+°° 3 
175.4/144 
0.04 



I 04+0.05 

4 - y4 -0.05 

II 61+ 113 

11.01_j 5g 

< 9 
8.54l° ; 2 
0.14_°;°? 
1 

79.5/72 
0.25 



r 

iVu 
3 [keV] 

7~z edge 

x> 



1.72: 



,+0.06 
a -0.06 

4.78l° : 8 7 « 
7.66±g:H 

o.4ot°:ii 

144.7/107 
8.9 x 10~ 3 



1 «s+o- 05 

l-88_o 05 
c\ 99+0.53 
Uiii -0.46 
+0.16 
' - d ' -0.16 

43+ 09 

u -^°_0.'09 

121.9/118 
0.38 



9 n7 +o.( 
z - u ' -0.( 



1.04 
J. 04 
r 4O+0.40 

7 oo+0. 18 
' • oo -0.17 

o.38l° o: °| 

133.9/108 
0.05 



9 19+O.O3 
^• ±z -0.03 

5 74+ - 32 

d - '^-0.31 

7-io±8:l| 
o.3ol° : ° 7 

186.4/146 
0.01 



1 QCC+0.06 
l.»O_0.06 

' 17 



11.2611H 
7 60+ ' 25 

' - DU -0.28 

21 + 09 

u - zl -0.09 

83.4/73 
0.19 



r 

iVu 
S zcdgc [keV] 

"7~zcdgc 

x> 



1 7 9+0.06 
1 -'^-0.05 



4.6: 



,+4.1 
-0.7 

<77 

7 68+ - 16 
' ■°°-o.ii 
+0.11 
0.10 



0.38: 
144.5/106 
7.8 x 10"" 



1 04+0.05 

15.69±g:I 
919+315 

7 qn+o-29 

I •OU_ 23 

S5+ 07 
u -°°-0.08 

112.8/117 
0.59 



9 1 o+0. 04 
z - ±o -0.05 
123 +3.94 

194j 



-4.11 

j+252 
-159 
7 1 7+O.I6 
'■- L '-0.09 
^+0.09 
-0.08 

126.4/107 
0.10 



0.36: 



2.16±°.8t 
10.37^° 2 

nq+209 
yo -76 

7 02 +0 11 

'• uz -0.12 

S1+ 07 
u - ol -0.07 

178.9/145 
0.03 



1 qc+0.08 
1 - 3J -0.07 



13.5 



+5.33 
3.66 
< 28 

7.08_°;?? 

o.i4tg : gg 

83.0/72 
0.18 



r 

iVu 
£a bs i[keV] 
o"absi [keV] 
EW absl [eV] 6 
£a bs2 [keV] 
o"abs2 [keV] 
EW abs2 [eV] 6 
X> 



1 74+0.06 
4 _ 89 +o.80 



8.05 
< 0.15 



0.73 

0.09 
0.07 



248±™ 

■0.19 
■0.18 
+0.17 
0.15 
104 
103 



9.78 
0.43 
461 
106.1/103 
0.40 



6.29: 
8.161 



1 8O+0.05 
L -° y -0.0B 
-,+0.55 
-0.52 
,+0.30 
-0.23 

o.3i±S-.lf 

380t 2 « 

W.94±l% 
9 nc+0.79 

Z.UO_o.53 

1008l 7Q -? 
121.8/114 
0.29 



9 n o+0.05 
z.uo_q 06 



5.41: 

8.281 



1+0.41 

-0.39 
j+0.41 
-0.38 

0.93_°-j° 
488l« 7 
12.04_J;?° 
2.95t 3 j? 
1 643+ 506 
121.0/104 
0.17 



9 1 1+0.04 
-0.03 
r 74 +0.32 
°- '^-0.31 

7 56+ ' 17 

'• OD -0.16 

S1+ 017 
U"jr_0.i4 

189± Io ° 
10.41_° : ?° 

9 nr+1.29 

864l 3 34 g 
179.0/142 
0.02 



I Q4+°-06 

II S5+ 1 ' 18 

__.OO_0 7 g 
o Q9+0-22 
°- oz -0.24 

0.48l° : 27 

+24 
25 
1.59 
0.73 



333 
13.51 



u.»__ . 92 
400+322 
^»»_320 

64.9/69 
0.62 



r 1 71+O.O5 1 nc+0.05 9 nq+0.07 9 1R+0-05 i Q9+O.O6 

1 1 -'*-0.06 ± - 3u -0.06 z - uo -0.06 z ' ±u -0.04 L -< i ^-0.06 

N H 4.71J:? 20.71?:? 5.80±g:°| 13.321IJ1 11.5±tt 

£ <77 498lg* <0.22 2281^ < 36 

£ absl [keV] 8.05!^ 7.76l°:l? 8.26l° ; 2 « 7.44« ; 1 10 8.39l°; 2 ^ 

a absl [keV] < 0.15 0.37l:°' 2 ? 0.95l°'* 3 0.27l:°' 1 1 2 0.37_o1? 

EWab.ifeV] 6 2411H 294ll 7 7 4 74l* 22 23ll°' 22 285l^ 7 

£a bs2 [keV] 9.78l°J? 10.681 J jl 12.07l 2 ' 3 ? 10.78±i"|| 13.86l ^ 

<7a bs2 [keV] 0.43l°i? 2.621^1 2.9ll?' 2 « 2.821?"^ l-67l ' 2 ^ 

EWab^feV] 6 455ll 8 7 7 963±||_ 15771^ 10581 387 476l 3 l 7 
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Table 2 — Continued 



Model" Parameter 


Values Epoch l d 


Values Epoch 2 d 


Values Epoch 3 d 


Values Epoch 4 d 


Values Epoch 5 d 




106.5/102 


106.4/113 


120.1/103 


163.3/141 


65.9/68 


p(x 2 M 


0.36 


0.66 


0.12 


0.10 


0.55 



a Model 1 is a power-law with Galactic absorption and intrinsic absorption (APL; XSPEC model wabs*zwabs*pow) ; 
Model 2 is a power-law with Galactic absorption, intrinsic absorption, and a notch absorber (APL+No; XSPEC model 
wabs*zwabs*notch*pow); Model 3 is a power- law with Galactic absorption, intrinsic absorption, and a notch absorber 
(IAPL+No; XSPEC model wabs*absori*notch*pow); Model 4 is a power- law with Galactic absorption, intrinsic absorp- 
tion, and an absorption edge (APL+Ed; XSPEC model wabs*zwabs*zedge*pow); Model 5 is a power-law with Galactic 
absorption, intrinsic ionized absorption, and an absorption edge (IAPL+Ed; XSPEC model wabs*absori*zedge*pow); 
Model 6 is a power-law with Galactic absorption, intrinsic absorption, and two absorption lines (APL+AL; XSPEC model 
wabs*zwabs*[pow+zgauss+zgauss]); Model 7 is a power-law with Galactic absorption, intrinsic ionized absorption, and 
two absorption lines (IAPL+2AL; XSPEC model wabs*absori*[pow+zgauss+zgauss]); All model fits include the Galactic 
absorption toward the source (Stark et al. 1992). TVh is the intrinsic column density expressed in units of 10 22 atoms cm~ 2 . 
All fits have been performed on the combined spectrum of images A, B and C of APM 08279+5255. 

b EW stands for equivalent width, which is defined as EW = f F " f Fe dE, where F c is the continuous flux and Fe is the 
flux in the absorber. 

C A11 absorption-line parameters are calculated for the rest frame. 

d All errors are for 90% confidence unless mentioned otherwise with all parameters taken to be of interest except absolute 
normalization. 

e P{\ 2 /v) is the probability of exceeding \ 2 for v degrees of freedom if the model is correct. 
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Table 3 

Results from spectral fits using XSTAR to the Chandra and XMM-Newton observations of APM 08279+5255. 
Model" Parameter 2 Values Epoch l d Values Epoch 2 d Values Epoch 3 d Values Epoch 4 d Values Epoch 5 d 



r 

A^Habsl 
log £absl 
-S'absl 

log £abs2 
2abs2 

r> 
P(x 2 M 



1 77 +0.04 
1 - 1 1 -0.04 

4 ^+0-5 

4 -°-0.3 

100 

4.otn 8 

3.2±0.10 

wtl 

q 9+0.2 
°- z -0.3 

2.4±0.1 
114.3/103 
0.21 



i 04+0.04 

1 - 9y: -0.04 

5 5+°' 7 

O.O_ 2 

17±2 

q 1+0.12 

3.0±0.1 

o c+0.7 
O.O_ 4 

2.0±0.1 
118.3/116 
0.42 



9 n7 +0.03 
^■ ul -0.02 
c q+0.2 
°-°-0.3 

39 
a fi +0.9 

3.1±0.1 

83 
o q+0.2 

1.7±0.1 
129.2/104 
0.05 



9 9n +0.03 
z ' z,u -0.03 



5.2^ 
12? 



,+0.2 
-0.4 
1+2.0 

q I +0.05 
— 0.17 

3.3±0.1 

q 1+0.13 

J - 1 -0.05 

2.3±0.1 
167.2/142 
0.07 



1 QQ+0.04 
1 - ao -0.04 

10.5±J"i 

8.5li; 8 5 

q q+0.7 

0.6 

2.94±0.10 

401^ 

o n+0.5 

1.48±0.07 
71.5/69 
0.40 



r 

-^direct 
A^Habsl 
log £absl 
-^absl 
-Kabsl 
AfHabs2 
log £abs2 
2abs2 
K a ,bs2 

X> 

P(x 2 M 



1 70+0-03 
L - ' y -0.04 

< lxl0~ 5 
100 

q Q + 0.07 

3-8-0.3 

3.2±0.1 

8l \ x 10~ 5 
22 

o q+0.4 

2.3±0.1 
5±1 x 10~ 5 
116.9/101 
0.13 



1.96 j 



H0.05 
-0.10 
c 7+0.3 
°-'-0.8 



2.9 



+ 1.0 

-2.2 



xlO" 



50 

q 1+0.14 
°- ± -0.07 
3.0 

7.8±\- 6 7 xlO" 8 
50l 7 2 
3.5 4 ~ 



- + 1.6 



(-0.6 

-"-0.6 
2.0 

xlO" 



116.2/114 
0.43 



2.0818-" 

c 9+0.20 
d -^-0.23 

1.2l^xl0- 5 
100 

q q + 0.06 

°-°-0.05 

3.04±0.1 

5.81°;? xHT 5 
100 
4 n+ 08 

4 - u -0.3 

1.6±0.1 
1.2tg-jxl(r 4 
116.2/102 
0.16 



2.32: 
5.5? 



>+o.i 

--0.09 
1+0.4 
'-0.3 

lAt° i xlO" 4 
7211° 

-,+0.07 
-0.07 



3.0! 
3.2 ±0.1 



1.2? 



;xio _ 



100 

q 1+0.3 
°- L -0.2 

2.3±0.1 
8I0.3XIO- 
156./140 
0.17 



1 QQ+0.03 
i - yo -0.04 

10.511° 
< 5xl0 -5 

is±!i 

q q+0.5 
o.o_ 07 

2.95 ±0.1 
l.ll^ xlO" 
100 

q O+0.2 

1.5±0.1 
6.3I0.3XIO- 
71.7/67 
0.33 



a Model 8 is a power-law with Galactic absorption, intrinsic neutral absorption, and two intrinsic ionized absorbers (APL ± 
2IA; XSPEC model wabs*zwabs*warmabs*warmabs(pow)). Model 9 is a power- law with Galactic absorption, intrinsic neu- 
tral absorption, and two partially covered intrinsic ionized absorbers (APL + PC*(2IA); XSPEC model wabs*zwabs*(pow 
+ warmabs*pow + warmabs*pow)). In the APL + PC*(2IA) model indirect, Knbsi and K a ,bs2 are the photon flux densities 
(photons keV -1 cm~ 2 s _1 ) at 1 keV of the one non-absorbed and two absorbed components, respectively. All model fits in- 
clude the Galactic absorption toward the source (Stark et al. 1992). All fits have been performed on the combined spectrum 
of images A, B and C of APM 08279+5255. Nh is the intrinsic column density expressed in units of 10 22 atoms cm~ 2 . We 
note that in several cases (especially for spectral fits using model 9) the value of the column density of the ionized absorber 
reaches the maximum allowable value in our XSTAR model of 1 x 10 24 atoms cm . The errors for several parameters are 
not listed because the fits in these cases could not provide any useful constraints. All listed errors are at the 68% confidence 
level. 
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Table 4 

Projected maximum outflow velocities, mass-outflow rates and efficiencies of outflows in APM 08279+5255 a . 



Epoch 


Instr. 


1>abal 

[c] 


log JVH(absl) 


M (absl) 


€K (absl) 


fabs2 

[c] 


log A r H (abs2) 


M (abs2) 
[Mo/i^Vr -1 ] 


eK (abs2) 


1 

2 
3 
4 
5 


ACIS S3 
EPIC pn 
EPIC pn 
EPIC pn 
ACIS S3 


o.isiroi 
.2o±r 3 4 

u - iz _0.02 
21 + 03 


22.9±0.4 
23.1±0.3 
23.2±0.3 
22.8±0.3 
23.0±0.7 


656t 8 4 f 

on i +1288 
oyj -_651 

209±»° 
609if 9 2 6 


o.outn 1 ! 4 
0.026±°;°£ 
O.04t° ;° 8 3 

n nna+0.004 
u - UUJ -0 .002 

n n3 +0 - 03 


n q« +0 - 02 

u " 3D -0.02 

46+ 08 

u '™-0.06 
n co + 0.07 

41+ 06 
61+ ' 07 

U.D1_ 04 


23.1±0.5 
23.4±0.3 
23.6±0.3 
23.4±0.4 
23.0±0.3 


1229 + 7 °f 

3141 + 899 

5508_ 3973 
2549+2593 

,, Q o + 1916 

ioao_ 1165 


u - 1 ' -0.11 

1 61 + 1 ' 89 
46+ - 47 

u -™-0.31 

62 +0 - 74 



a The estimated values of the outflow properties were based on fits that assumed an absorbed power-law model with two absorption lines. The 
values of M and ek are obtained by equations 1 and 2, respectively, assuming Mbh ~ 10 12 /i.£ ^0 C see §3-2) and Lw — 7 X 10 15 /j^ Lq (Irwin 
ct al. 1998; Reichers ct al. 2009). 
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Table 5 

The minimum and maximum energies and velocities of the high-energy absorption features in APM 08279+5255. 



Epoch 


Instrument 


-Em in 


-Em ax 


^min 


V max 






[keV] 


[kcV] 


[c] 


[c] 


1 


ACIS S3 


8.05±0.17 


10.64±0.27 


0.18+0.02 


0.43+0.02 


2 


EPIC pn 


7.54±0.74 


15. 04+1. 22 


0.12+0.09 


0.67+0.04 


3 


EPIC pn 


6.42±0.82 


17.94+2.96 


0.04+0.04 


0.76+0.07 


4 


EPIC pn 


6.94±0.35 


16.31+1.17 


0.04+0.04 


0.71+0.04 


5 


ACIS S3 


7.36±0.56 


15.35+2.12 


0.09+0.07 


0.68+0.07 



a i5 m in and -E ma x are estimated from spectral fits using model 6 of Table 2. 
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Fig. 1. — Ax residuals between the best-fit Galactic absorption and power-law model and the Chandra, ACIS spectra of 
APM 08279+5255. This model is fit to events with energies lying within the ranges 4.5-10 keV. The arrows indicate the best-fit 
energies of the absorption lines of the first and second outflow components for epoch 1 (top panel) and epoch 5 (lower panel) 
obtained in fits that used model 6 of Table 2. 
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Fig. 2. — Ax residuals between the best-fit Galactic absorption and power-law model and the XMM-Newton pn spectra of 
APM 08279+5255. This model is fit to events with energies lying within the range 4.5-10 keV. The arrows indicate the best-fit 
energies of the absorption line of the first outflow component for epoch 3 (middle panel) and epoch 4 (bottom panel) obtained in 
fits that used model 6 of Table 2. 
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Fig. 3. — 68%, 90%, and 99% confidence contours between the normalizations of the absorption lines at i? a b s i and i? a bs2 and 
the respective energies -Eabsi and i? a bs2 for the Chandra 2002 (top panel) and Chandra 2008 (bottom panel) observations of 
APM 08279+5255 assuming fit 6 of Table 2. The 99% confidence contour of component abs2 for epoch 5 is not well constrained 
and the 95% contour is shown instead. 
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Fig. 4. — 68%, 90% and 99% confidence contours between best-fit energies i? a bsi and _E a bs2 and energy widths <j a b s i and <7 a bs2 
of the absorption lines for the Chandra 2002 (Epoch 1) and Chandra 2008 (Epoch 5) observations of APM 08279+5255 assuming 
fit 6 of Table 2. The 99% confidence contour of component ahs2 for epoch 5 is not well constrained and not shown in this plot. 
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Fig. 5. — 68%, 90% and 95% x 2 confidence contours of _E a bs versus cr a b s of the first absorption line in epoch 3 (dotted line) and 
epoch 4 (solid line) assuming the APL + 2 AL model (see model 6 of Table 2). The 95% confidence contour for epoch 3 is not 
well constrained and not shown in this plot. 
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Rest-Frame Energy (keV) 
2.445 4.91 9.82 24.55 49.1 




Observed-Frame Energy (keV) 

Fig. 6.— The spectra and best-fit models (model 6 of Table 2) of the XMM-Newton observations of APM 08279+5255. We note 
the significant variability of the X-ray BALs (extending between w 1.5-3.6 keV) and continuum between epochs 3 and 4 that are 
separated by only ~ 3.3 days (proper-time). 



22 



CHARTAS ET AL. 




Fig. 7. — The spectra and best-fit models (model 6 of Table 2) of the Chandra observations of APM 08279+5255. We note the 
significant variability of the X-ray BALs extending between » 1.5-3.6 keV. 
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Fig. 8. — The ratios of the (a) observed and (b) best-fit-modeled (using model 6 of Table 2) flux densities Fi/Fj between the 
three XMM-Newton observations of APM 08279+5255 listed in Table 1, where i and j represent the epochs compared, At is the 
observed-frame time between epochs and F is the flux density in units of counts s _1 keV -1 . Note that a significant change of flux 
density is observed between all observations over the entire spectrum of APM 08279+5255. The spectral change between epochs 
3 and 4 is particularly remarkable since it occurs over a period of only ~ 3.3 days (proper-time) implying a size-scale of ~ 7.4 x 
10 15 cm which is comparable to nsco = 4.5 x 10 15 cm for the case of a Schwarzschild black hole in APM 08279+5255. 
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Fig. 9. — The ratios of the (a) observed and (b) best-fit-modeled (using model 6 of Table 2) flux densities F5/F1 between the two 
Chandra observations of APM 08279+5255 listed in Table 1, where Fi and F5 are the flux densities for epochs 1 and 5 respectively. 
At is the observed-frame time between epochs. Note that a significant change of flux density is observed near the region of the Fe 
BALs (~ 1.5-3.6 keV). 
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Fig. 10. — Maximum outflow velocity, i>max, observed from APM 08279+5255 as a function of the X-ray photon index, T. « ma x and F were 
derived from fits to the spectra of APM 08279+5255 with a model that included an absorbed power-law and two Gaussian absorption lines 
(model 6 of Table 2). Errors shown are at the 68% confidence level. Data shown with circles, a star, and squares are obtained from Chandra, 
SuzaJcu and XMM-Newton observations respectively. 
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Fig. 11. — The bound-free, Mq, and bound-bound, Ml, components of the force multiplier are shown as a function of the ionization 
parameter. Force multipliers are calculated for SEDs with photon indices of T = 2.1 (solid lines) and V = 1.7 (dotted lines). In the left 
panel we have assumed no absorbing shield, whereas, in the right panel the soft and hard SEDs have been attenuated by an absorber with 
log Nn = 22.8. 
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Fig. 12. — Maximum outflow velocity as a function of the angle between the outflow direction and our line of sight. Shown are the outflow 
velocities for the most likely resonance transitions responsible for the observed X-ray BALs in APM 08279+5255. The outflow velocities were 
calculated using the relativistic Doppler formula for the case of the maximum rest-frame energy of the X-ray BALs of E max = 17.9 keV (see 
Table 5) observed in epoch 3. 



